Quantum mechanical and information theoretic view on classical glass transitions 
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Using the mapping of the Fokker-Planck description of classical stochastic dynamics onto a quantum Hamil- 
tonian, we argue that a dynamical glass transition in the former must have a precise definition in terms of a 
quantum phase transition in the latter. At the dynamical level, the transition corresponds to a collapse of the 
excitation spectrum at a critical point. At the static level, the transition affects the ground state wavefunction: 
while in some cases it could be picked up by the expectation value of a local operator, in others the order may 
be non-local, and impossible to be determined with any local probe. Here we instead propose to use concepts 
from quantum information theory that are not centered around local order parameters, such as fidelity and en- 
tanglement measures. We show that for systems derived from the mapping of classical stochastic dynamics, 
singularities in the fidelity susceptibility translate directly into singularities in the heat capacity of the classical 
system. In classical glassy systems with an extensive number of metastable states, we find that the prefactor of 
the area law term in the entanglement entropy jumps across the transition. We also discuss how entanglement 
measures can be used to detect a growing correlation length that diverges at the transition. Finally, we illustrate 
how static order can be hidden in systems with a macroscopically large number of degenerate equilibrium states 
by constructing a three dimensional lattice gauge model with only short-range interactions but with a finite 
temperature continuous phase transition into a massively degenerate phase. 



I. INTRODUCTION 

Order and disorder are fundamental concepts in under- 
standing phases of matter. A classification of states of mat- 
ter is possible according to the symmetries that are broken in 
the condensed state? 1 which can be detected through the non- 
vanishing expectation value of a local observable: the order 
parameter. Crystalline order for instance is characterized by 
the breaking of full translational symmetry into a regular pe- 
riodic lattice. Glasses, however, do not have a (spatially) local 
order parameter that can be defined from a single snapshot 
of the microscopic configuration of the physical degrees of 
freedom, say, the position of the constituent atoms. [There 
are order parameters in replica type theories, but they require 
correlations across replicas. 2 ] Hence, glass order does not fall 
within the Landau paradigm of classification via symmetry 
and local order parameters. 

It has been argued for many years that the glass transition 
might have a purely dynamical origin, and that it is not accom- 
panied by a thermodynamic (static) transition. This scenario 
prompts then a fundamental question: Is it possible to define 
accurately a dynamic phase transition, and to understand why 
static order appears to be absent? 

In this paper, we show that a positive answer to this question 
can be found by exploiting the mapping between the stochas- 
tic dynamics of a classical system and a quantum mechan- 
ical system at zero temperature. The mapping is based on 
a well-known symmetrization of the classical transition ma- 
trix^, provided it satisfies detailed balance, and on the pro- 
motion of the classical configurations to orthonormal basis 
states of a Hilbert space. This yields a correspondence be- 
tween the classical system and a quantum Hamiltonian in the 
preferred configuration basis, with two fundamental proper- 
ties: (i) all quantum ground state correlators of diagonal op- 
erators in the preferred basis are equivalent to classical corre- 
lation functions; and (ii) the energy spectrum of the Hamilto- 



nian equals the spectrum of relaxation rates in the stochastic 
classical system. In particular, the longest equilibration time 
scale in the classical system is inversely proportional to the 
splitting between the ground state and the first excited state 
of the quantum Hamiltonian. (See Ref. |4j for a constructive 
formulation of this correspondence.) In essence, the map- 
ping merges a d-dimensional classical system with stochas- 
tic dynamics (i.e., with an additional 'time dimension') into a 
d-dimensional quantum system. In the quantum mechanical 
language, concepts in quantum phase transitions and quantum 
information can be used to formulate precisely what a dynam- 
ical transition is in the original classical system, and provide 
new avenues to study these transitions. 

In conventional approaches to dynamical transitions, one 
often relies on time dependent concepts, such as metastable 
states around free energy minima. These concepts require the 
ability to distinguish between short and long relaxation time 
scales, that characterise for instance the relaxation of the sys- 
tem within a basin, and transitions between basins, respec- 
tively. This in turn introduces an intrinsic difficulty in com- 
paring the behaviour of finite size systems with respect to the 
thermodynamic limit. 

Note however that the quantum Hamiltonian contains all 
the information about the original classical stochastic process, 
and, as such, one is not required to look at time dependent 
quantities. In recent years, a view of stochastic dynamics as "a 
kind of thermodynamics in space-time"— was developed using 
large deviation functions^ In this theory, one weighs trajec- 
tories by a "Boltzmann measure" that couples to their 'activ- 
ity', thus favouring slow vs fast relaxing modes in the system, 
depending on the sign of the coupling constant. In the quan- 
tum mechanical language, the weights amount to perturbing 
terms in the Hamiltonian which break detailed balance, and 
transitions are studied as a function of the coupling constant. 
In a sense, large deviation functions can be seen as a way to 
probe the physics of the quantum system by looking at how it 
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responds to external perturbations. Relations between large 
deviations in classical systems and quantum critical points 
have been discussed recently in Ref. [4g, and further elabo- 
rated in Ref. S3 

In this paper, we take a different point of view. By commit- 
ting to the study of the quantum mechanical system per se, we 
show that one can detect, characterise and possibly classify 
glass transitions as peculiar static quantum phase transitions. 

If the classical system undergoes a dynamical transition, 
the gap must close in the associated quantum model, leading 
to a degeneracy that persists throughout the incipient phase. 
Such gapped-to-gapless transition signals therefore a (ther- 
modynamic) quantum phase transition that provides the pre- 
cise definition of a dynamical transition in the classical sys- 
tem. (Generically, soft modes can be present also in the fast- 
relaxing, high-temperature phase of the system. In Sec.HITlwe 
discuss how these soft modes ought to be treated separately 
from the spectral collapse responsible for a glass transition.) 
The appearance of a non-vanishing Edwards-Anderson order 
parameter— across the transition translates directly into a di- 
vergent static local susceptibility in the quantum system. 

One of the great conceptual difficulties in defining a dy- 
namical transition is the fact that there seems to be no local 
order parameter differentiating the phases on the two sides of 
the transition. A glass looks locally and instantaneously like a 
liquid, much as a spin glass looks like a paramagnet. Here we 
argue that using the quantum mechanical mapping one can do 
away with local order parameters in the classical configuration 
basis, and one can effectively investigate static signatures of 
a dynamical glass transition using off-diagonal operators and 
non-local measures, which are more natural in the quantum- 
equivalent formulation of the stochastic problem. 

A quantum phase transition is signaled by changes in the 
ground state wavefunction. For example, if there is a local 
order parameter, the transition can be picked up by the ex- 
pectation value of local operators. These operators need not 
be diagonal in the preferred basis given by the classical con- 
figurations. For example, in a spin system where the classi- 
cal basis is the z-direction, ordering in the x-direction in the 
quantum system amounts to a 'dynamical' order in its clas- 
sical counterpart. The quantum mechanical language gives 
immediate access to operators written in a form that is not 
diagonal in the preferred basis, and therefore it broadens the 
field of search for local order parameters that might capture a 
glass transition. 

Moreover, quantum phase transitions can be studied us- 
ing measures that do not rely on the existence (and there- 
fore, on any a priori knowledge) of an order parameter. For 
instance, concepts from quantum information theory, such 
as fidelity susceptibility^^ and entanglement entropy^ have 
proven useful in detecting and characterising transitions into 
exotic phases with non-local order parameters (e.g., topologi- 
cal^ ordered phasesi^ili). Motivated by the difficulty in find- 
ing any static order accompanying dynamical transitions, we 
propose to employ such measures in this context. 

The fidelity is constructed from the overlap of the wave- 
function at two infinitesimally close values of the coupling 
constant that drives the transition^ It can be expressed 



straightforwardly in the basis where the classical configura- 
tions are defined. Here, we show that, for quantum Hamilto- 
nians derived from the stochastic dynamics of a classical sys- 
tem, the quantum fidelity susceptibility is proportional to the 
heat capacity of the classical system. A singularity in the fi- 
delity susceptibility at the phase transition translates therefore 
into a singularity in the heat capacity of the classical system. 
The latter is a feature that is indeed present in several theoreti- 
cal and experimental examples of glass transitions, 15 but upon 
which there is so far no consensus in the literature. 

Another tool to study quantum phase transitions, even ex- 
otic ones without a local order parameter, is the von Neumann 
entanglement entropy. Scaling properties of the von Neumann 
entropy as a function of the size of the subsystem, for instance, 
can reveal a correlation length that diverges at the phase tran- 
sition in topologically ordered systems. It is thus interesting 
to examine what entanglement entails for the quantum system 
derived by mapping the dynamics of a glass. Here we discuss 
how the von Neumann entanglement entropy of a bipartition 
of the systems is sensitive to the properties of the two classi- 
cal subsystems, and to the way they influence one another. For 
instance, in the case of glassy phases with an extensive num- 
ber of non-relaxing modes (in the thermodynamic limit), we 
show that the entanglement entropy obeys the (conventional) 
area law both above and below the transition, but with a pref- 
actor that jumps across the phase transition. As we discuss, 
one can also use it to detect a diverging correlation length at 
the dynamical transition. 

The paper is organised as follows. We begin by re-deriving 
typical (time dependent) quantities used to characterise glass 
transitions in the quantum mechanical language. For exam- 
ple, the Edwards-Anderson order parameter^ defined as the 
appropriate limit of a connected autocorrelation function, can 
be straightforwardly written in quantum mechanical language 
to show that a phase with a non-vanishing ggA is in fact char- 
acterized by a divergent static susceptibility. The susceptibil- 
ity is computed with respect to a perturbation to the mapped 
quantum Hamiltonian that cannot be written as an equivalent 
perturbation of the original classical system. This perturba- 
tion is "quantum" in nature, akin to the one used in large devi- 
ation theory. [The familiar reader might want to move quickly 
through Sec. |il]and Sec.imil 

We then enter the discussion of static quantum mechani- 
cal measures that can be used to investigate glass transitions. 
In particular, we focus on quantities that probe phase transi- 
tions without resorting to local order parameters: the quan- 
tum information theoretic concepts of fidelity susceptibility 
(Sec. lIVI ) and entanglement entropy (SeclVT). 

Finally, in Sec. [VI] we construct an explicit example where 
we can analytically show that there is a static phase transi- 
tion accompanying a dynamical one, in a system where the 
order is hidden. This example, a generalization of the three- 
dimensional gonihedric spin model, has a finite temperature 
transition into a phase with a macroscopic number of equi- 
librium states. If endowed with gaugeable disorder, it is im- 
possible to select a single minimum-energy configuration by 
applying a local field. 

In order to ensure a clear flow through the central results, 
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several parts of the paper have been moved to the appendix. In 
Appendix[A] we discuss the nature of the collapsing states at 
the transition. In Appendix iBl as an example of the mapping 
of classical stochastic dynamics to quantum Hamiltonians, 
we carry out explicitly the construction for the Sherrington- 
Kirkpatrick model, and we explain the connection between 
the Parisi equilibrium minima and the degenerate quantum 
states. In AppendixICl we show how one can construct varia- 
tional wavefunctions for the case of symmetry-broken states, 
and how ergodicity breaking can be related to their vanishing 
variational energies. In AppendixlDl by writing the gonihedric 
model in terms of dual variables, we construct an example of a 
(toric-code-likei£) quantum system with a macroscopic topo- 
logical ground state degeneracy. 



II. REVIEW OF THE MAPPING OF STOCHASTIC 
PROCESSES TO QUANTUM HAMILTONIANS 

The fact that a transition matrix W in a Markov process 
(obeying detailed balance) can be made into a symmetric ma- 
trix by means of a similarity transformation has long been ex- 
ploited to either find the eigenvalues and eigenvectors of Wr 
or to gain insight on quantum Hamiltonians by turning them 
into classical dynamical systemsJi 

In Ref. @ it was shown that this correspondence is indeed 
one-to-one onto a class of quantum Hamiltonian representa- 
tions, dubbed Stochastic Matrix Form (SMF) decompositions, 
that essentially encompasses all quantum systems devoid of 
the sign problem. 

In this mapping, the temperature T = 1//3 of the classi- 
cal system enters the quantum Hamiltonian as a coupling con- 
stant. In particular, equilibrium correlators in the classical sys- 
tem at temperature T map onto equal time ground state (GS) 
correlators in the quantum system, and the Markov relaxation 
rates become the spectrum of the quantum Hamiltonian. 

The study of eigenvalue degeneracy in dynamical systems, 
suggested by Kac amongst others, has long been exploited in 
trying to characterise the appearance of long-lived metastable 
states, and glassiness (see for instance Refs. 18,191 and refer- 
ences therein). 

Consider a system whose configurations are labeled by C, 
and where the probability Pq (t) of the system being in con- 
figuration C at time t is governed by the Markov equation 



dt 



P c (t) = Y,Wc,c>Pc(t), 



(2.1) 



where the transition matrix Wc.c satisfies probability con- 
servation Wc,c = — 2c#C^ C ' C '' an d detailed balance 
%,e"^/ T '= W a ,ce- E c>/ T . 

It is convenient to introduce the following vector notation. 
We can consider Pc (t) as being the real component of a vec- 
tor \P(t)) in the span of the orthonormal basis {\C) , VC}, 
P c (t) = (C\P(t)). Eq. <|2J) becomes then 



dt 



\P(t)) = W\P(t)), 



(2.2) 



which admits solutions in the form 

\P(t))=exp(wt) |P(0)>, 



(2.3) 



where we used the fact that W does not depend explicitly on t. 
Note that Wc,c is rea l but not necessarily symmetric, hence 
W needs not be Hermitian, and one has to distinguish between 
left and right eigenvectors. 

The right eigenvectors are given by 



W 



(2.4) 



where Eq — 0, e n > for n ^ (as ensured by Perron- 
Frobenius theorem for any irreducible finite system size), and 



cm 



(2.5) 



[For convenience, in the following we shall assume that the re- 
laxation rates are labelled in ascending order, £o = < £i < 
£2 < • ■ •] We can then decompose \P(t)) into normal modes, 

and use Eq. (H2J to obtain |P(t)) = E„ ««(0) e~ £ "* 

In this language, we can express the ensemble average of any 
given observable O as 

(0(t)) th = Y,O c P C (t) (2.6) 

c 

= J2(C\d\P(t)) = (E\d\P(t)), (2.7) 
c 

where we defined |E) = J2c \ C )> and O \ C ) = °c 10- Note 
that the state (S| = ^ c (C\ when bracketed with any prob- 
ability state gives (H\P(t)) = 1, i.e., the probability is nor- 
malised at all times, ensured by the fact that (S| = (ip^ 

is a left eigenvector with zero eigenvalue. Decomposing into 
normal modes, we obtain 



(2.8) 



If £1 remains finite in the thermodynamic limit, for t — > 00 
only the slowest mode £0 = survives and one arrives at the 
thermodynamic equilibrium expression 



(0(^oo)) th = (E|0 



-E C /T 

~z ' 



We are interested in understanding in a precise way how the 
spectral properties of the operator W control the relaxation 
properties of the system. To this end, we investigate the two- 
time autocorrelation function: 

C(t + r,t) = {0(t + r)0(t)) th 

= Y,°C' p c^C'(r) O c Pc(t) 

c,c 



^(c'\de WT \c) (c\6\p{t) 



c,c 



Y,\Oe WT 0\P(t)) 



(2.9) 
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where Pc^c ( T ) is the conditional probability that the system 
be in configuration C at time t, given that it was in configu- 
ration C at time t = 0. 



■li- 



on 



of W are not neces- 



Note that the right eigenvalues 
sarily orthonormal. Thanks to the detailed balance condition, 
it is possible to construct a real symmetric matrix Hsmf that 
has the same eigenvalues as W by means of a similarity trans- 
formation using the matrix Sc.c = cxp (~/3Ec/2) Sc.c'- 



Hsmf — —S 1 WS, 



(2.10) 



with (-ffsMF)cc = (Hsmf)c,C following from detailed bal- 
ance. 

Labeling the eigenvectors of H by we take advantage 
of the fact that (n\m) = S njn and J2 n \ n ) ( n \ = ' to express 
Eq. (fZ9]i as 

C(t + r,t) = e Wr 6\P{t)) 

= ( s l OSe-" SUF T \n) (n| S^O \P(t)) 
= (0| S-'OS |n) (n| S^O \P{t)) , 

n 

(2.11) 

where we used the conservation of probability condition, 
namely that (E| W = 0, to rewrite (S| = (0| -S' -1 . If we 
take the limit f ^ oo, 

C(t) = lim C(t + t, t) 

t—>oo 

= 53 e" £ " T (0| S^OS \n) (n\ S^O \P(oo)) 

n 

= J3 e" £ " T (0| S^OS \n) (n\ S^OS |0) , (2.12) 

n 

where we used the fact that |P(oo)) = S\0)M Classical 
thermodynamic observables are by definition measured in the 
preferred basis of classical configurations {|C)}. Therefore 
the operator O is diagonal in this basis, it commutes with S 
(S- 1 OS = 6), and 



III. EDWARDS-ANDERSON PARAMETER 

If we order the eigenvalues £o — < e\ < £2 < • ■ ., 
we can write an upper bound of the connected autocorrelation 
function as 

<e- ElT J2(0\6\n)(n\6\0) (3.1) 

n 

= e~ eiT (0| 6 2 |0) . (3.2) 

Unless £1 — > in the thermodynamic limit N — > 00, the 
connected correlation function decays to zero exponentially 
fast in time (ri = 1/ei), independently of the choice of the 
(bounded) observable O. Therefore, any dynamical transition 
leading to 'long' (i.e., non exponential) decay must be accom- 
panied by a vanishing spectral gap in the SMF Hamiltonian. 

The appearance of vanishing relaxation rates in stochas- 
tic processes can be related to several different factors: for 
instance, critical slowing down at symmetry breaking phase 
transitions, Goldstone modes when a continuous symmetry is 
broken, or diffusive modes in paramagnetic phases with con- 
served quantities. Quite generally however these phenomena 
are considered distinct from glassiness in that they lead to 
parametrically large relaxation time scales (i.e., the gap closes 
as a power of system size). Glassiness on the contrary is char- 
acterised by a gap that vanishes exponentially in system size. 
In this paper we focus on the latter phenomenon. Without 
loss of generality and for ease of discussion, we neglect slow 
modes that vanish as a power law of system size - the reader 
can consider for example the case of a system with discrete 
degrees of freedom and no conserved quantities, illustrated 
qualitatively in Fig. Q] [Compare with a similar assumption 



A Energies £ n 



c(r)=xy 



O|0) 



(2.13) 




Notice that the r.h.s. in the above equation is precisely the GS 
quantum imaginary-time autocorrelation function (7 quant (t) 
obtained from the operator O.— 

The connected part of the correlation function can then be 
written as 



C c (r) = (0(t + T)0(t)) th -(Oy th 

2 



71^0 

C c quant (r). 



>n\O\0) 



(2.14) 
(2.15) 



FIG. 1 . Qualitative illustration of the spectral collapse in the associ- 
ated quantum system which is expected to occur at a glass transition 
in the original classical system with discrete degrees of freedom and 
no conserved quantities. 

in Ref. [TH] As a result, the high-temperature paramagnetic 
phase in the systems discussed below always relaxes expo- 
nentially fast in the thermodynamic limit, and the associated 
quantum Hamiltonian is in a gapped phase. 

Dynamical glassy phases are typically characterized by a 
finite value of the Edwards-Anderson order parameter, ^ea = 



5 



lim T _ ) . 00 C c (t)i& In quantum mechanical notation. 



Qea(0)= |ho 

n£T>,n^0 



(3.3) 



where T> is the set of eigenvalues e n that vanish in the limit of 
N — > oo (i.e., they collapse onto the GS). 

Let us postpone for the moment further considerations on 
how a non-vanishing qea arises in a quantum system, and let 
us look at its physical significance. As we pointed out earlier, 
the classical autocorrelation function in Eq. ( 12.131 l is equiv- 
alent to the imaginary-time quantum autocorrelation func- 
tion C quant (r) ee (0\e HsMFT de- HsMFT 6\0). By tempo- 
rally integrating the connected part of the quantum correlator, 
Q uant (r) ee C"» uant (r) - | (0| 6 |0) | 2 , one obtains the local 
static (zero-frequency) susceptibility of the quantum system 
at zero-temperature, 



x loc (^ = 0) 



dr C quant (r) 



E 

n#0 



n|O|0 



(3.4) 

The Edwards-Anderson order parameter <7ea, defined in 
Eq. ( 13.31 ). corresponds to the long time limit r —> oo of the 
connected correlator C quant (r). Therefore, the appearance of 
a non-vanishing value of (?ea coincides with the onset of a 
divergent local static susceptibility M 

Notice that the susceptibility in Eq. (13.4b measures the re- 
sponse of the quantum SMF Hamiltonian to a perturbation that 
couples to the operator O = J^c \C)Oc(C\. The value of ^ea 
is thus related to the second order correction in the GS en- 
ergy of the perturbed system, Hsmf + AO. We point out here 
that this perturbation is quantum mechanical in nature, and 
in general it is not equivalent to some classical perturbation 
Eq — > Eq + SEc in the original system. 

In conclusion, a dynamical transition in the original classi- 
cal system becomes a well-defined thermodynamic quantum 
phase transition in the associated Hamiltonian ffgMF- This 
is a transition from a gapped to a gapless phase, detected by 
a non-vanishing order parameter (qea), and accompanied by 
the appearance of a divergent static susceptibility. In other 
words, a dynamical phase transition can be precisely identi- 
fied as an equilibrium phase transition in the quantum system. 

Conceptually the quantum perturbation that gives rise to the 
local static susceptibility discussed in this section is equiva- 
lent to those considered in large deviation functions^ Indeed, 
in the latter one perturbs the transition matrix W that controls 
the dynamics of the classical process, W — > W + sK. The 
perturbation K is intended to couple differently to slow and 
fast relaxing modes in the system. By tuning the coupling 
constant s away from s = 0, one thus favours/disfavours the 
low-lying eigenstates of W. Upon performing the similar- 
ity transformation in Sec. [II] K maps onto a static perturba- 
tion K = —S KS to the quantum mechanical Hamiltonian 
H = — S~ 1 WS. So long as K can be expressed as a local 
term in space and time in its classical formulation, the quan- 
tum operator K is local, and one can qualitatively understand 



the behaviour of the system as a function of the coupling con- 
stant s, 7 as follows. At high temperature the quantum sys- 
tem is gapped (see above for a discussion of this assumption). 
Therefore, varying s can induce a transition in the system only 
at finite values of s (if at all present). On the other end, be- 
low a glass transition, there is a large number of degenerate 
lowest-lying states. At least in finite size systems, the degen- 
eracy is not exact, and the operator K is likely to split the 
degeneracy, inducing a first order transition at s = 0. In or- 
der to make this argument rigorous, the thermodynamic limit 
ought to be properly accounted for, which is beyond the scope 
of the present paper. (Note that this behaviour occurs also at 
classical symmetry breaking transitions.) 



IV. FIDELITY SUSCEPTIBILITY 

The information theoretic concept of fidelity has been re- 
cently used to detect quantum phase transitions without a pri- 
ori knowledge of any order paramateri^iiS It has also proven 
useful to investigate the nature of a critical point, allowing 
one to derive scaling exponents and other universal proper- 
ties. This technique provides a general approach to study a 
dynamical phase transition, that encompasses the case where 
no local order parameter is found even at the quantum level 
(e.g., topologically ordered systems^). 

The fidelity overlap is defined as 

T(/3, 6/3) ee (M(3 - 5/3/2) | ifotf + 6/3/2)) , (4.1) 
from which one obtains the so-called fidelity susceptibility 



lim 

<5/3— s-0 



sp 2 



(4.2) 



One of the features of the classical to quantum correspon- 
dence is that the GS wavefunction of the quantum system is 
known exactly for all values of (3, which allows for a direct 
calculation of xf{P)- F° r a generic SMF Hamiltonian, the 
GS wavefunction can be written as 



IM/3)>=E 



exp(-/3£ c /2) 



\C), 



(4.3) 



where Z{J3) = ^2 C exp(— /3Ec/2) is the partition function 
of the classical model. From it, we obtain that the fidelity 
susceptibility \T °f m e SMF quantum system is proportional 
to the heat capacity Cy of the original classical system, 



xAP) 



2 \nZ((3) 



Ad/3 



(4.4) 



where we used the fact that E = — d\n Z((3)/d(3, and that 
Cy = dE/dT is the heat capacity of the classical system. 

Singularities in the heat capacity (not necessarily a diver- 
gence) are therefore in one-to-one correspondence with sin- 
gularities in the GS fidelity susceptibility. Whilst there is no 
rigorous proof that all quantum phase transitions in systems 
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with local Hamiltonians give rise to singularities in \jr, no 
counter examples are known, including quantum disordered 
and topologically ordered systems.— This is suggestive that a 
generic dynamical phase transition in a local classical system 
is detected by the fidelity susceptibility as a thermodynamic 
transition in the associated quantum system, and it is likewise 
signalled by a necessary singularity in the heat capacity of the 
original classical system. 



guage. Nonetheless, as we discuss below, it captures the prop- 
erties of metastable states in the original system, which can 
only be defined by a distinction between "short" and "long" 
time scales in the classical time-dependent language. 

For concreteness, consider a standard bipartition of the de- 
grees of freedom of the system into two subsets: subsystem 
A, a bubble of finite radius R, and subsystem B, as shown 
in Fig. |2] Using the density matrix constructed from the GS 



V. ENTANGLEMENT ENTROPY TO PROBE GLASS 
TRANSITIONS AND DIVERGING LENGTH SCALES 

The von Neumann entanglement entropy has recently been 
applied quite extensively as a tool to study quantum systems. 1 1 
It provides an unbiased measure, in the sense that it does not 
focus on a particular order parameter or local operator, to 
probe properties that are not accessible via standard correla- 
tion functions among a fixed number of degrees of freedom. 
Because it is not hinged on order parameters, the entangle- 
ment entropy can be a means to uncover "hidden" order that 
is not easily detectable otherwise. For example, it has been 
used to detect topological orders in systems where no local 
order parameters existi22r— 

One measure of entanglement is block entanglement, ob- 
tained by bipartitioning the system into two blocks A and B, 
and computing the von Neumann entropy Sab for the reduced 
density matrix after tracing out the degrees of freedom in B 
(or A). When the system is short range correlated, the en- 
tanglement entropy Sab obeys the so-called area law,— i.e., 
Sab oc l AB ' wriere £ab is the length scale of the bound- 
ary between the subsystems and d is the dimension of space. 
At criticality, it is possible to have logarithmic corrections to 
the area law, as is the case in one dimension (ID).— How- 
ever, this is not always the case in higher dimensions, and it 
is possible to have no corrections to the area law even when 
the correlation length diverges^ A correction to the area law 
near a critical point is not the only signature of a quantum 
phase transition. In the case of a topological phase transition, 
subleading corrections to the entanglement entropy that cap- 
ture the topology of the bipartitions do show a clear change as 
the system crosses transitions between topological and non- 
topological phases. 12 Moreover, one can use the scaling of the 
(topological) entanglement entropy as a function of the size 
of the bipartitions to define a growing correlation length that 
diverges at the trasition, even if the system is devoid of a local 
order parameter. 

Here we propose that the entanglement entropy for the 
mapped quantum system can be used as a tool to study glass 
transitions that are hard to probe otherwise, say, using equal- 
time correlation functions. Although the entanglement en- 
tropy is generally difficult to compute, the point that we make 
in this section is that the behavior of the entanglement entropy 
can provide a formal way to uncover glass transitions without 
a priori knowledge of any order parameter, and without the 
need of ad hoc definitions of time-scale-dependent free energy 
basins. We remark that computing the entanglement entropy 
amounts to a static calculation in the quantum mechanical lan- 




FIG. 2. Bipartition of a system into subsystems A and B. 



wavefunction Eq. ( 14.31 ). 



m = \MP))(MP)\ 



HP) 



-p(E c +E c ,)/2 



\C){C'\ 



(5.1) 



and tracing over subsystem B to get the reduced density ma- 
trix Pa{P) = Trs/5(/3), one obtains the entanglement en- 
tropy S A b{P) = -Tta[pa(P)^Pa{P)}- The entaglement 
entropy for a pure state, like the one in Eq. (14.3b . is symmet- 
ric, i.e., one arrives at the same result by computing Sab = 
-Tr B [p B {p)lnp B (P)]. 

The computation of Sab from Eq. d5.lt was carried out 
explicitly in Sec. V of Ref.l31[ 

Sab(P) = PF A + /3F Q - /3F AuB + (E d ) th (5.2) 



where /3F A = -lnZ£, f3F B = -InZg, and [3F AuB = 
— In Z. The partition functions Zp\ and Z B are computed 
with Dirichlet boundary conditions on the interface between 
the two subsystems, i.e. they are computed by pinning the 
configurations at the boundary of the bipartition and running 
the sums over configurations in the bulk of A or B, respec- 
tively. That is to be contrasted with Zp and Zp, where 
the partition function is computed with free boundary con- 
ditions, simply summing over all configurations in A or B 
without boundary constraints. Notice that the total partition 
Z{f3) = Zp UB is naturally computed with free boundary con- 
ditions, as there is no boundary in the total system A U B. Fi- 
nally, (E d ) th is the average energy of the interface between 
A and B (i.e., the sum of all terms in Eq that involve simul- 
taneously degrees of freedom in A and in B) computed with 
the measure P(C) = exp(-/3E c )/Z(p). 
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It is instructive to rewrite the boundary entropy in the fol- 
lowing form, 

S A b(P) = P(F A - Ff) + f3(F B - FE) 

+ (3Fl + /3F£ - (3F AuB + (3 (E d ) th . (5.3) 

We can then define the free-boundary-condition part of the 
entanglement entropy 



where F 9 = F AuB — F A — F B = /T 1 ln(exp (f3E 9 )) 
and we can introduce the notation 



(5.4) 



ur 



AF A}B = F A B - F A B 



A,B 



A,B 



so that 



Sab(P) = P {AF A + AF B ) + S% B (0). 



(5.5) 



(5.6) 



The term S AB (/3) can identified with the cumulant- 
generating function for the fluctuations of the boundary en- 
ergy, 



iJ 



AB 



(j3) = In (e 



j(E 8 -{E B ) tb ) 



th 



(5.7) 



The second cumulant is a measure of the heat capacity for 
the boundary degrees of freedom. S AB (j3) generically obeys 
the area law away from quantum critical points. If there is a 
singularity in the heat capacity of the system, then this term 
becomes singular and picks up the phase transition, much like 
the fidelity susceptibility in SeclIVI 

However, as we argue in the following, the entanglement 
entropy S A b in Eq. ( 15.61 ) is able to detect glass transitions 
even in subtle cases when there are truly no singularities in the 
specific heat / quantum fidelity. This can be seen by consid- 
ering the other terms in Sab, namely AFa,b, which contain 
information on the effects of the boundary conditions on the 
free energy of subsystems A and B. 

Let us adopt a description of the phase space of subsys- 
tem A in terms of metastable free energy basins, commonly 
used in the literature of glassy phenomena. The goal here 
is to argue that the von Neumann entropy exhibits signatures 
of glassy phenomena, and it can be used to determine grow- 
ing length scales and other characteristic features. On the 
other hand, however, such time-scale-dependent description 
in terms of metastable states is not needed to compute / study 
the von Neumann entropy in the quantum mechanical lan- 
guage, and we argue that Sab provides a static measure to 
investigate those same glassy phenomena, without the need to 
justify a metastable state description of the system. 

Consider the bipartition (A, B) of the system illustrated in 
Fig. [2] Let us first allow the whole system to equilibrate, 
and then let us freeze subsystem B and let A thermalize with 
fixed boundary conditions. In the high-temperature liquid 
phase the fixed boundary conditions play a marginal role, and 
the free energy of A is characterised by a single minimum, 
which is not significantly different from the one that obtains 
for free boundary conditions. On the other hand, glassiness is 



typically accompanied by the appearance of distinct free en- 
ergy minima, and the partition function of A can thus be de- 
composed accordingly into separate contributions. Following 
Bouchaud and Biroli^, we make the coarse assumption that 
all these contributions have the same bulk free energy, and dif- 
fer only by a boundary energy term. Moreover, we take this 
boundary term to be the same for all contributions (Eg) but 
for one that is favoured (Eg — AE) by an amount AE > 0. 
If we label these minima in the free energy of A by indices 
ai, i = 1, . . . ,Af, with interfacial energy E a .^ 1 ~ Eg, and 
E ai ~ Eg — AE, then under the phenomenological assump- 
tions introduced above it is possible to write 



1 7 A 



i=l 



1 



Af 

CAT — 1) + eP AB 



Af 

1 + e ^E-s- 



(5.8) 



where S* = InM is the configurational entropy of A. 

From Eq. ( 15.81 ). we can obtain the behavior of the contribu- 
tion AFa to the von Neumann entanglement entropy of the 
associated quantum system across the glass transition temper- 
ature. Let us focus on the case where the configurational en- 
tropy S* becomes extensive in the glassy phase. Above the 
glass transition, the liquid phase (where Af = 1, and S* = 0) 
quickly relaxes to minimize the energy strain due to the fixed 
boundary conditions, and AFa ~ E ai = Eg — AE. Below 
the transition, the appearance of long relaxation time scales 
prevents the system from moving between the Af minima. For 
large enough R, the extensive entropy S* ~ R d dominates 
the exponential in Eq. d5.8t , and AFa = Eg. 

This behaviour becomes evident if we take A to be half of 
the entire system, and we consider the thermodynamic limit. 
In this case, AFa = AF B . Irrespective of S* being fi- 
nite or extensive, the area law is obeyed (AFa oc i? d_1 ), 
but with different coefficients (namely, (Eg — AE)/R d ^ 1 vs 
Eg/R d ^ 1 ). Therefore, at least within the scenario where the 
number of minima scales extensively with (sub)system size, 
the glass transition can be captured by a sudden change in 
slope of the von Neumann entropy of a bipartition of the sys- 
tem as a function of the interface area, as illustrated in Fig. [3] 
(Note that this conclusion is reached within a mean field pic- 
ture that neglects possible additional logarithmic corrections 
to the area law at the transition temperature.) 

We stress that, although we used a metastable state descrip- 
tion to obtain the behaviour of AFa, the entanglement en- 
tropy is a static measure, that is defined in terms of the quan- 
tum ground state of the mapped Hamiltonian, i.e, it is related 
to the classical equilibrium state. Therefore, the entangle- 
ment entropy can be used independently of whether we have a 
proper understanding of the metastable states in the classical 
system. 

In addition to its behavior across the transition, the depen- 
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A SabM 





FIG. 3. Qualitative illustration of the behaviour of Sab /A in the 
thermodynamic limit as a function of T, if we take, say, subsystem 
A to be half of the entire system. 



FIG. 4. Qualitative illustration of the behaviour of AFa as a func- 
tion of area A(R), as predicted in the scenario discussed in Ref.l32l 
below the glass temperature T g . 



dence of the von Neumann entanglement entropy on the size 
of the bipartition can be used to detect a growing correlation 
length £ that diverges at the transition. Here, as in Ref. [lH £ 
is identified as the crossover size of A such that S* ~ (3AE 
(recall that S* ~ R d and AE ~ R*- 1 ). The value of £ de- 
pends on the temperature both explicitly (j3) and through the 
(intensive) configurational entropy s* = S*/R d . Within the 
metastable state picture, it follows from Eqs. ( 15. 6t and d5.8t 
that for R < £ one finds AFa ~ E* d - AE = E ai (i.e., the 
energy difference dominates over the entropic gain), while for 
R 3> £ one finds AFa ~ -Eg (i.e., the entropy dominates, and 
the free energy difference between the fixed and free bound- 
ary conditions is given by the boundary energy of the majority 
of the configurations of A). Thus, a plot of AFa (and there- 
fore of Sab) as a function of the interface area A(R) ~ 
exhibits a kink around R ~ £, where the slope changes from 
F Ql /^ to£^/.4 (see Fig. 0). 

We would like to stress the following important point. The 
length scale defined by Bouchaud and Biroli required a care- 
ful definition of the metastable states, based on the existence 
of a regime of time scales sufficiently large compared with 
the equilibration time within each metastable state, but small 
compared to transitions between the states. Again, here we 
used their results just to argue that the entanglement entropy 
does capture this growing length scale, and exhibits a dis- 
tinctive behaviour across the transition. However, the entan- 
glement entropy per se does not require any discussion of 
metastable states and time scales. The entanglement entropy 
computed from Eq. (15. 6t is a static measure of the equilibrium 
state of the system. Our point is that this static but non-local 
measure can be used to probe/reveal a "hidden" static transi- 
tion, and the associated growing correlation length. 

In summary, we propose that the entanglement entropy can 
be used as a way to probe glass phases and phase transitions 
without bias towards order parameters. Within the scenario 
used for instance in Ref. [3^ we showed that the entangle- 



ment entropy obeys the area law, with a prefactor that changes 
abruptly across the glass transition. It can also be used to iden- 
tify static growing length scales that in the classical language 
require a time-dependent formulation in terms of configura- 
tional entropiesr^i^ or point-to-set correlation functions^ 
The entanglement entropy that we propose to use is a com- 
pletely static measure that can serve as a probe of non-local 
order, and thus it can be a tool to reveal an underlying static 
transition associated with the dynamical glass transition. 



VI. HIDDEN ORDER IN AN IDEAL GLASS PLAQUETTE 
MODEL 

In this section we discuss a class of models with short-range 
interactions that exhibit a large number of degenerate low- 
est energy configurations, and we consider the nature of the 
energy barriers between such configurations at low tempera- 
tures. In particular, for one of the models we show that the 
existence of a thermodynamic phase transition can be argued 
by mapping it to essentially decoupled 2D Ising layers. Be- 
low the phase transition, the energy barriers between the many 
energy minima grow exponentially with system size, and the 
existence of an equal number of slow relaxing modes can be 
shown rigorously by means of a variational calculation involv- 
ing the associated SMF Hamiltonian. 



A. Generalised gonihedric model 

Consider Ising spin variables o\ defined on the sites i = 
{ix, iyi iz) of a cubic LxLxL lattice (with periodic boundary 
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conditions imposed), with energy^! 

E = —J X y <7i Ci+x <7j + x + y 0"i+y 

i 

~Jyz CT i ^i+y ^i+y+z °"i+z 
i 

— Jzx ^2 &\ fi+z fi+z+i fi+x ) (6-1) 
i 

where x, y and z are the unit vectors of the cubic lattice. The 
system contains interactions among four spins at the vertices 
of square plaquettes, with coupling constants that depend on 
whether a plaquette is parallel to the xy, yz, or zx plane. We 
focus on the case where 



J - J' 



j 



(6.2) 
(63) 



which has an anisotropy between the coupling constants for 
vertical and horizontal plaquettes (with respect to the xy 
plane). 

Two particular limiting cases of this model are equivalent 
to systems already studied in the literature. For J = 0, J' ^ 
the system is equivalent to L decoupled copies of the 2D 
square plaquette model.^^i This model is known to exhibit 
an activated behaviour, with time scales growing as the ex- 
ponential of the inverse temperature^- 6 The model does not 
exhibit a dynamical transition, and the longest relaxation time 
scale diverges only in the zero temperature limit. 

The isotropic limit J = J' corresponds to the gonihedric 
model.^2, This model has been argued to undergo a first order 
thermodynamic transition, using numerical simulations and 
cluster mean field arguments^ Moreover, numerical simula- 
tions point to the existence of a finite temperature glass transi- 
tion in the neighbourhood of the thermodynamic transition.— 
However, no analytical approach has been successful at con- 
firming the presence of such dynamical transition thus far. 

The case we shall consider here is the limit J ^ 0, J' = 0, 
which we refer to as the anisotropic gonihedric model. As 
we show below, this model has a thermodynamic transition at 
T c = 2 J j log(l + -\/2), which is the value of the 2D Ising tran- 
sition temperature for a nearest-neighbor interaction J (see 
Refs. |5o| and H3 for a study of the system based on transfer 
matrix considerations). In contrast with a single Ising plane 
that has two degenerate minima, this system has a number of 
degenerate energy minima that scales as 2 i_1 . 



1. J' = - kinetically constrained model with a finite 
temperature transition 



Let us write the anisotropic gonihedric model as a gauge 
theory, introducing Ising degrees of freedom 9n$\ defined on 
the links (or bonds) (ij) between nearest neighbor sites i and 
j of the cubic lattice (see Figure 0, 



^{ij> = (TiCTj. 



(6.4) 



The new Ising degrees of freedom are subject to the hard con- 
straint that the product of four 9 spins along the edges of any 



given plaquette must be equal to 1 (gauge constraint). Under 
periodic boundary conditions, there are additional constraints 
on the 9 spins, since products along lines winding around the 
system (say, parallel to the x, y, or z axis) must also equal 1. 



-I 



P 



<ij> 



FIG. 5. Illustration of the mapping from a (solid circles) to 8 (solid 
squares) spins, for a plaquette p of the cubic lattice. 

Because of the gauge constraint, one can actually work 
solely with 9 variables on the vertical bonds between sites i 



and i + z, with centers at b = (i x , i 



i z + 



1/2): 



7b — O-b-z/2 Cb+z/2- 



(6.5) 



There is a 1-to-l mapping between the original a spin con- 
figurations, and the (constrained) configurations of 6*b spins, 
provided that we specify the values of the o\ spins on a given 
reference xy plane, e.g., i = (i x ,i y ,0). In the limit of interest, 
J / and J' = 0, the energy of the system can be expressed 
purely in terms of vertical 6*b spins, 



E 



b 



b+x 



7b Cb4 



(6.6) 



The plane of a spins required by the mapping behaves as a 
separate paramagnetic contribution to the system, which can 
be disregarded. 

In the 9 spin language, the limit J' = of the generalised 
gonihedric model in Eq. d6.lt can be recognised as a collec- 
tion of decoupled 2D Ising models in disguise [see Eq. ( 16.61 )1. 
This is not a trivial result, and in the original a spin formula- 
tion the J' = model, while being clearly anisotropic, is far 
from being factorisable into decoupled 2D layers. 

Note that the decoupling is not exact. Indeed, the prod- 
uct of theta spins along vertical lines Jlf =i ^b=(i x ,i H ,i z +i/2) 
must equal one because of the nature of the mapping (gauge 
condition). However, we do not expect the parity constraint to 
affect the properties of the system in the thermodynamic limit, 
as it simply reduces the number of independent 9 degrees of 
freedom by L 2 . 

Below the 2D Ising transition temperature T c = 
2 J j log(l + y/2) the system orders, and each layer (which we 
can label by k = 1, . . . , L) acquires an (independent) expec- 
tation value M k = (% x , lH ,fc+i/2)> = ±M(T/T C ) (apart from 
the parity constraint, enforcing ]l£=i ^=(^,^,^+1/2) = 
+1). There are therefore 2 L_1 minima that correspond to 
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all relative magnetizations of the layers. This is an exam- 
ple of a system without disorder but with a true thermody- 
namic transition into a phase with many minima, whose num- 
ber scales subextensively with system size (namely, exponen- 
tially in L = ^N). 

In a 2D Ising ferromagnet below T c , the system settles in 
one of the two minima (spontaneous symmetry breaking), and 
the energy barrier between them grows with the linear size of 
the system. As a result, the time scale for the system to mi- 
grate from one minimum to the other grows exponentially in 
L (see for instance AppendixICt. The same is true for the sys- 
tem in question, except that the number of minima separated 
by energy barriers ~ L scales as 2 L . We therefore expect 
to observe an equal number of independent slow relaxation 
modes in the transition matrix of the system (assuming lo- 
cal a spin dynamics), with rates ~ cxp(— L), which vanish 
in the thermodynamic limit. In Sec. lVIBl we explicitly show 
the existence of these slow modes by means of a variational 
approach in the associated SMF quantum language. 

Whereas the 9 language allows one to reduce the system 
to decoupled 2D Ising planes, simple dynamical processes in 
the a spins translate into coordinated rearrangements in the 
9 spins - a feature that is characteristic of kinetically con- 
strained models. In the present case, local dynamics in the a 
spins are bound to couple the dynamics of different 9 planes, 
reflecting the fact that the system is in truth three dimensional. 
For instance, single a spin flip events translate into flipping 
two vertically adjacent 9 spins, that is, spins that belong to 
two adjacent planes: 



(Ti — > —0\ 




1-8/2 
i+z/2 



(6.7) 



This is an example of a kinetically constrained model with a 
finite temperature thermodynamic transition. The coordinated 
dynamics can affect the coarsening taking place within each 
2D layer of 9 spins. For example, in order for a domain wall 
in one of the layers to move, 9 spins must be flipped either 
in the layer immediately above or below, as a consequence 
of Eq. (16.7b . If these spins happen to be within an ordered 
domain, as opposed to a boundary between domains, they in- 
duce an energy cost to the domain wall motion in the adjacent 
layer, and therefore coarsening becomes an activated process. 
We can summarize the features of this model as follows: 

(I) A a spin can be written as a product of vertical 9 spins and 
a a in the reference plane. The average of the product scales as 
[±M(T/T c )] d , where M(T/T C ) e (0, 1) is the magnitude of 
the average magnetisation of a spin, and d is the distance to the 
reference plane. Therefore, (a) —> exponentially fast away 
from the reference plane. The system is a spin liquid, with no 
long-range order in the cr-spins. [See also Ref. [34], where it 
is shown using transfer matrix arguments that the correlation 
length of the system remains finite across the transition.] 

(II) Nevertheless, {w^) = (0 i+2/2 ) = ±M(T/T c ),\f i z = 
1, . . . , L, and there are sub-extensively many minima (2 L_1 ). 
Provided that the dynamics in the a spins are local, the system 
is characterised by an equal number of slow relaxation modes, 
whose rates vanish exponentially ~ cxp(— L) in the thermo- 
dynamic limit. 



(Ill) In addition to these exponentially large time scales, 
the two-time average (<7i(i) cr i+2 (i) <ri(t w ) cr i+2 (t w )) (= 
(@i+z/2(t) @i+z/2{tw)}) shows at least t/t w scaling (coarsen- 
ing within the planes), if not an even slower aging behaviour 
due to the coordinated inter-plane dynamics. 



2. J' 7^ - the relation to the kinetically constrained square 
plaquette model, and confinement 

Suppose that one starts from the limit J' ^ 0, and J = 0, 
i.e., truly decoupled 2D square plaquette layers. As was 
shown in Ref. [3a it is well understood that the dynamics be- 
come slow as temperature is reduced. Single a spin flip up- 
dates change the sign of all four adjacent plaquettes in the 
plane, so in order to move an isolated defective plaquette (one 
where the product of the four a spins at the corners is nega- 
tive) requires the creation of extra pairs of defective plaque- 
ttes, each at a cost J' . However, pairs of defective plaque- 
ttes can move freely across the system, and defects can there- 
fore be annihilated or separated by processes that encompass 
only finite energy barriers, repeated for a number of steps that 
scales linearly with the distance between the defects. As a re- 
sult, the system exhibits a conventional activated behaviour, 
and the time scales diverge only in the zero temperature limit. 
In other words, defects in the 2D square plaquette model are 
deconfined, as only a finite energy is required to separate them 
infinitely far apart. 

Let us now discuss the effect of turning on a small J. While 
the physics within each xy plane remains the same, flipping a 
a spin generates 8 further defective vertical plaquettes. Pro- 
cesses that efficiently separated, say, pairs of horizontal defec- 
tive plaquettes in the J = case, now acquire an additional 
energy cost that scales with the number of flipped spins. In- 
deed, moving defects within a plane leaves behind a wake of 
flipped vertical plaquettes. This "string" connecting the de- 
fects has an energy cost proportional to J times the separation, 
and the defects become linearly confined. A confinement- 
deconfinement transition occurs when the temperature is large 
enough to overcome the confining energy scale J, and the 
physics of the square plaquette model is recovered. Therefore, 
the model with J ^ is qualitatively different from J = 0. 

If we start with J ^ and decrease the ratio J'/ J, we ex- 
pect the same confined physics to survive at sufficiently low 
temperatures (depending on the value of J). What could per- 
haps change is the nature of the confinement-deconfinement 
transition at some critical temperature T c ( J, J') (indeed, as 
a function of Jj J', we observed using numerical simula- 
tions what appears to be a change in the transition character 
from second to first order). Therefore, we conjecture that the 
low temperature phase of the anisotropic gonihedric model 
(J' = 0) discussed in the previous section remains essentially 
unchanged when the ratio J'/ J is increased from zero. In 
support of this conjecture, we verified numerically that the lo- 
cation of the (putative) first order transition in the isotropic 
J' = J gonihedric model in Ref. [37] is in good agreement 
with the expected transition temperature that is obtained from 
2D Ising considerations, = 3T c /2 = 3J/log(l + V2) 
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(see also Refs. I41p4l Note that J = 1/2 in those references, 
as well as in Ref.l37l). 



A (gaugeable) random version of the generalized gonihedric 
model 



Here we discuss yet another generalization of the model in 
Sec. IVI A II That model has an extensive number of min- 
ima, namely 2 L_1 for a system of linear size L. The or- 
der parameters were (cicrj+s) = (#i+z/2) = ±M(T/T C ), 
Vi z = 1, . . . , L (there was one constraint, so only L — 1 
of these are independent). Suppose that in Eq. (16. U we take 
J xy = and generalize the other couplings so that 



E 



i 

_ J\x a i a i+z °"i+z+x O-i+x 
i 



with ,P ZX = ± J and J* z = ± J random, but "gaugeable", sat- 
isfying the condition J\ x J y + X J*+ y J yz = J 4 . In this non- 
frustrated case the signs of J^ x z ^ 2 and Jy Z z ^ 2 in Eq. ( 16.8b 
can all be gauged out, and the thermodynamic behaviour 
of the system is the same as in Sec. IVI A II undergoing a 
phase transition precisely at the 2D Ising transition temper- 
ature T c = 2J/log(l + \/2), with a 2 L_1 degenerate "or- 
dered" phase. However, in this model it is impossible to write 
down the order parameter in terms of local products of a vari- 
ables, without using non-local strings of products of the ran- 



dom variables J^ x z ^ 2 



± J and J,, 



example of a gaugeable glass 



42,43 



b-z/2 



± J. This is an 



B. SMF Hamiltonian and variational excited states 

Let us now consider the J' = limit of the generalised 
gonihedric model discussed in Sec. IVI A 1 1 from the quan- 
tum mechanical perspective (illustrated in greater detail for 
the Sherrington-Kirkpatrick model in App. |B}. [For similar 
considerations on the isotropic gonihedric model, we refer the 
reader to Appendix IDll 

Recall that the Ising degrees of freedom <Tj live on the sites 
i of a cubic lattice (with periodic boundary conditions, for 
simplicity), and the energy is proportional to the sum over all 
vertical plaquette products, that is 



E — — J ^ (9b $b+x + Ob Ob+y) , 



(6.9) 



where we introduced the Ising variables # b=i+ £/ 2 = °~i 
as in Eq. d6.5t . 

Assuming single spin flip Glauber dynamics, and follow- 
ing the steps outlined in App. [B] (see Ref.|4J), we arrive at the 



associated quantum SMF Hamiltonian 



-HsMF=y^ 



1 



i 2 cosh 



2 



i 



b 2 cosh 
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AE U (6.10) 



and the corresponding GS wavefunction is given by the super- 
position 



exp 



{¥E b (*£' 



b+x 



'b "b+y 



>)} 



(6.11) 

where Z = £ c cx P {/?J ^ 9 h+i + 9 h 6 b+ -)}. 

As discussed in Sec. IVI A 11 the behaviour of the ground 
state of this model is effectively described by decoupled 2D 
classical Ising layers, which is particularly evident in the 9 
spin language. As a result, there is a critical value T c where 
the gap closes in the quantum system, and the fidelity suscep- 
tibility exhibits a singular behavior (namely, ~ Cl slng , 
see Sec. HVI i. diverging logarithmically as the transition is ap- 
proached. 

Below T c , the classical system becomes massively degen- 
erate, with energy barriers between the lowest energy states 
that scale with the linear size of the system (see also Ref.H). 
These degenerate states give rise to an equal number of slow 
relaxing modes, and therefore to an equal number of low-lying 
excited states in the associated quantum system, whose energy 
tends to zero in the thermodynamic limit. 

Using the variational approach explained in detail in Ap- 
pendix [C] we can explicitly find an upper bound to a num- 
ber of low-lying eigenstates of Hsmf for T < T c equal to 
the number of thermodynamic energy minima in the original 
classical system. The upper bound tends to zero in the ther- 
modynamic limit, thereby proving that the dynamical classical 
system has many relaxing modes whose decay rates vanish in 
the thermodynamic limit (These are equivalent to the degen- 
erate eigenstates discussed in Ref. [HI where the equality be- 
tween the numbers of free energy minima and of low-lying 
dynamical states is proven rigorously. See also Ref. Qjl for a 
characterisation of metastable states in finite size systems in 
the original classical language.) 

The order parameter distinguishing the different free en- 
ergy minima below T c is the vector (Mi, . . . , Ml), where 
Me = ±\M(T/T C )\ is the average magnetisation in each 
plane. [Recall that a non-vanishing magnetisation of the 9 
spins along one plane corresponds to a non-vanishing expecta- 
tion value of the sum of products of nearest-neighbour a spin 
pairs, one immediately above and one immediately below that 
9 spin plane (see Fig.[6j)]. Following the steps in AppendixICl 
we expect low-lying variational excited states of the form 
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Notice that the argument of the hyperbolic tangent in C{ n } (C) 
is the z component of the magnetisation of the 8 spins in 
a given plane, Mt = Y^ L m , S z . /on-. Since 

n£ = 0, 1, for ^ = 1, . . . , L, and given the constraint that 
Yle=i must be positive, there are 2 L_1 independent varia- 
tional wavef unctions. 

By construction, (?/>{„} |-0o) = 0, and (ip{n}\tp{ n '}) = 0, 
V {n} 7^ {n'}, and eigenstates of H^mf must exist with en- 
ergies A {n} < (ip{ n }\H SMF \ip {n} ). Given that the thermo- 
dynamic properties of the anisotropic gonihedric model are 
controlled by the classical 2D Ising model, one can straight- 
forwardly follow the steps in Appendixlcl to find that an upper 
bound to A{„j is given by 



{n} (C)\ 2 - 0* {n} (C)0 {n} (C u 



ih 



(Pi 2 : 



'0(L 3 



< 



L 



T>T C 

(6.13) 



where the configuration Ck is obtained from the configuration 
C upon changing the sign of the 9 spins at sites k and k + z; 
(. . .) th is a classical thermal equilibrium average (used here 
in the mathematical sense, i.e., (Mi) th = even if the sys- 
tem undergoes spontaneous symmetry breaking); and a is a 
positive constant. 



VII. CONCLUSIONS 

In this paper we use a well-known mapping of classi- 
cal stochastic processes onto quantum Hamiltonians^^ (see 
Ref. for a detailed constructive approach) to argue that dy- 
namical glass transitions can be interpreted in the quantum 



mechanical language as static zero-temperature phase transi- 
tions where a number (exponentially large in system size) of 
excited states collapse onto the ground state. The quantum 
mechanical perspective allows one to accurately define what 
a dynamical glass transition means, and could provide new 
avenues to understand the consequences of the transition on 
the static properties of the system. Here we proposed to use 
concepts from quantum information, such as entanglement en- 
tropy and fidelity susceptibility, as tools to help uncover "hid- 
den" order in glasses, which may not be accessible through 
local order parameters. 

In mapping the stochastic dynamics of classical systems to 
quantum Hamiltonians, the relaxation spectrum of the classi- 
cal system corresponds to the excitation spectrum of the quan- 
tum model. Any interesting (i.e., non-exponentially decaying) 
dynamical behaviour of the former requires necessarily the 
vanishing of one or more relaxation rates. That is, the appear- 
ance of modes, other than the equilibrium distribution, which 
fail to relax during the stochastic process. Understanding such 
dynamical features in the classical system is tantamount to un- 
derstanding the behaviour of the spectrum in the associated 
zero-temperature, static quantum system4£ 

In the past few years a great deal of knowledge has been 
accumulated in the study of quantum phase transitions. In 
particular, it has become clear that there are transitions for 
which it is not possible to define local order parameters. Nev- 
ertheless, there are ways of detecting and characterizing such 
transitions without resorting to an order parameter. One re- 
cent tool is the concept of fidelity, which is constructed from 
the overlap of the ground state wavefunction at two infinites- 
imally close values of the coupling constant that drives the 
transition. Here we showed that the fidelity susceptibility, 
in quantum systems derived via the mapping from classical 
Forker-Planck evolutions, is precisely the heat capacity of the 
classical system. Therefore, we argue that, if the fidelity picks 
up quantum transitions in local Hamiltonians, then a dynami- 
cal transition in a classical system with local energies and lo- 
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cal dynamics must be accompanied by a thermodynamic sig- 
nature (singularity) in the heat capacity. We note that there are 
examples of disordered systems where dynamic and thermo- 
dynamic transitions differ, such as the fully connected p-spin 
glass modelr^ however, these are non-local systems. While 
there is no rigorous proof that all quantum phase transitions 
give rise to singularities in the fidelity, we know of no counter 
example in local Hamiltonians, including quantum disordered 
and topologically ordered systems. Indeed one can view our 
results as forcing a marriage between quantum information 
and non-equilibrium glassy physics, with strong implications 
from one field into the other. One is thus pressed to accept 
that either A) there exist strange undetectable quantum phase 
transitions in local Hamiltonians without local order parame- 
ters or even fidelity singularities, or B) there is a true thermo- 
dynamic glass transition, without necessarily any local order 
parameters, but with a heat capacity singularity nonetheless. 

Another tool that has proven effective to probe exotic states 
without a priori knowledge of an order parameter is the von 
Neumann entanglement entropy. Here we studied its behavior 
in quantum states constructed via the mapping from stochastic 
dynamics. We showed that, in the case of glassy phases with 
an extensive number of non-relaxing modes (in the thermo- 
dynamic limit), the entanglement entropy of a bipartition of 
the system obeys the (conventional) area law both above and 
below the transition, but with a prefactor that jumps across 
the phase transition. In addition, scaling properties of the von 
Neumann entropy as a function of the size of the subsystems 
can be used to reveal a correlation length that diverges at the 
dynamical phase transition. 

Finally, in this paper we introduced an example of a non- 
disordered spin model that exemplifies well the concepts aris- 
ing from the classical to quantum mapping, and it exhibits 
several of the features typical of a "glassy" spectral collapse. 
The model, discussed in Sec. [VI] is an example of a sys- 
tem with purely local interactions that can undergo a thermo- 
dynamic transition at finite temperature into a phase with a 
(sub)extensive number of equilibrium states. The model has 
2 L ~ 1 degenerate equilibrium states, and its (gaugeable) dis- 
ordered version is an exactly solvable ideal spin glass, with 
purely local interactions and a (sub)extensive number of sta- 
tistically equivalent (from a local point of view) equilibrium 
states. 

Understanding exotic zero-temperature quantum phase 
transitions is not obviously simpler than studying glass tran- 
sitions. The mapping discussed in this paper is not a magic 
wand, but rather a change in perspective. Yet sometimes a 
change in perspective is what one needs to gain new insight 
on a long standing problem. 
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Appendix A: Nature of the collapsing states 

As we discussed earlier, any dynamical transition requires 
a collapse of relaxation rates, i.e., of eigenenergies in the as- 
sociated quantum Hamiltonian H^mf- Let us focus here on 
the nature of these collapsing states, in particular their re- 
lation to spontaneous symmetry breaking, broken ergodicity, 
and glassiness. 

Consider first the conventional case of a classical thermo- 
dynamic symmetry breaking transition with a local order pa- 
rameter. Below the transition temperature, the classical free 
energy landscape develops distinct minima separated by en- 
ergy barriers ~ L a , for some exponent a < d. The number 
AT of such minima is controlled by the broken symmetry, and 
it is typically finite (i.e., it does not grow with the size of the 
system): if the broken symmetry is discrete, there are a fi- 
nite number of such states; if the symmetry is continuous, the 
manifold of degenerate states is labelled by a continuous vari- 
able. This in turn leads to JV low lying eigenstates of H^mf, 
with energies that scale as E\ ~ e~ aL (recall that the ground 
state energy is by construction Eq = 0). Take for instance the 
Ising ferromagnet in Appendix ICl 

The nature of the spectrum above these exponentially de- 
generate states depends on the dynamics governing the relax- 
ation within a minimum. For instance, coarsening leads to a 
power law e cx ~ L~ z level spacing. Notice that one can argue 
for Goldstone modes (e ox ~ L~ b ) in the continuous symmetry 
breaking case, but one can still have algebraic (in 1 /L) energy 
levels even in discrete systems, as the Hsmf spectrum de- 
pends on the dynamics (consider for instance diffusive modes 
in presence of conserved quantities). 

The exponentially large times required to switch between 
broken symmetry minima give rise to broken ergodicity in the 
thermodynamic limit. However, the presence of such diverg- 
ing time scales is immaterial, because they can be observed 
only in a fine tuned system. The presence of a local order 
parameter means that an infinitesimal local perturbation^ to 
the energy of the classical system produces a finite separation 
between the energies of the different minima - with conse- 
quent removal of the diverging time scales. In studying the 
behaviour of the system, we are thus justified to pick one of 
the symmetry-broken states over the others. 

Such a relation between symmetry breaking order parame- 
ters and broken ergodicity is well known. There are however a 
number of other scenarios that lead to a collapse of relaxation 
rates in classical Markov processes. 

(i) Systems with broken symmetry where the degener- 
acy scales with system size - This is the case for the example 
without quenched disorder discussed in Sec. [VI] An important 
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difference between this case and the one where the degener- 
acy does not scale with system size is that it is not possible 
to split all the degeneracy by applying an infinitesimal local 
perturbation. However, it is still usually possible to select a 
unique GS, and one could argue that the physics in the end is 
not much different from a conventional spontaneous symme- 
try breaking transition. To illustrate this, let us consider the 
example in Sec. I VI A II where the system behaves as decou- 
pled 2D Ising ferromagnetic layers. An infinitesimal uniform 
magnetic field h lowers the energy of the uniformly magne- 
tised state by an amount > hL 2 with respect to all other states 
that are degenerate for h = 0. Note that, crucially, hL 2 di- 
verges in the thermodynamic limit for any infinitesimal but 
fixed value of h. 

(ii) Disordered glassy systems - If we consider well 
known examples of glassy systems such as the Sherrington- 
Kirkpatrick model, discussed in App. [B] a degeneracy 
emerges at low temperatures, at least if we adopt Parisi's pic- 
ture, which scales with system size. The minima are statisti- 
cally equivalent from a local perspective, and unlike the case 
of magnetic ordering in an Ising ferromagnet, there is no local 
operator that we can apply uniformly across the system to lift 
the degeneracy. Clearly if one knew the lowest energy spin 
configuration in one of the minima, a magnetic field could 
be tailored to favour that precise configuration energetically. 
Whilst this is indeed a perturbation given by the sum of local 
operators, the values of the field are random but fixed specifi- 
cally for each and every spin in the lattice. Such random field 
contains highly non-local information, and we feel that this 
should not be called a local perturbation. 48 

(ii) Systems with gaugable (i.e., not frustrating) disor- 
der - The leading difference between glass transitions and 
symmetry breaking transitions - namely, the ability to select a 
unique GS by means of infinitesimal local perturbations - can 
be removed by introducing disorder in a gaugeable fashion 
(see for example the model in Sec. lVI A~3l l. A gaugeable dis- 
order does not spoil the thermodynamic behaviour, but with- 
out full knowledge of the gauge transformation, one can no 
longer distinguish the degenerate states below T c using local 
perturbations. For a model such as the anisotropic gonihedric 
one, with a GS degeneracy that scales with the size of the sys- 
tem, the addition of gaugeable disorder yields a model that 
becomes indeed very similar to a disordered glassy system. 
One might then wonder whether the gaugeability of the dis- 
order does in fact constitute a material difference between the 
models, or if the physics is ultimately the same from a quali- 
tative point of view. 

(iv) Systems with topological degeneracy - For complete- 
ness, it is interesting to compare the scenarios above (where 
the collapsing spectrum was drawn from classical stochastic 
processes) with yet a different case drawn from quantum me- 
chanics, although there is no a priori connection with clas- 
sical stochastic processes in this case. Consider a zero tem- 
perature quantum system which undergoes a phase transition 
into a topologically ordered phase. In this phase, the system 
develops a topological degeneracy which is ultimately non- 
local. Not only does any infinitesimal local perturbation fail 
to produce a finite lifting of the degeneracy, but actually any 



local operator must have identical expectation values across 
all the GSs, and vanishing matrix elements between any two 
degenerate states. For example, the reader can think of modi- 
fications of the toric code Hamiltonian J£ As shown in Ref.[l2i 
one can indeed construct SMF-like quantum Hamiltonians 
that exhibit topological order in regions of their phase dia- 
grams. In systems where the number of topologically degen- 
erate states scales with the size of the system (e.g., in models 
similar to those in Sec. IVII briefly outlined in Appendix |D), 
the spectral collapse at a transition into the topological phase 
is qualitatively that of a perfect glass; so perfect that it yields 
no measurable consequences! Indeed, the classical free en- 
ergy minima are exactly identical from a statistical point of 
view, and the diverging time scales for the system to go from 
one minimum to another are immaterial to all locally measur- 
able quantities (e.g., the Edwards-Anderson order parameter 
vanishes on both sides of the transition). In order to see a di- 
verging time scale, one would need, say, to compute non-local 
(loop-loop) autocorrelation functions. 

To conclude, when the collapse of the energy spectrum is 
due to the spontaneous symmetry breaking, or to the emer- 
gence of topological order, the system goes through a quan- 
tum phase transition in ffgMF into a phase with a manifold of 
exponentially degenerate states AE ~ e~ aL . In this paper, 
we argue that also a dynamical (glass) transition becomes a 
well-defined static phase transition in the associated quantum 
mechanical language. 

Appendix B: The Sherrington-Kirkpartick model, revisited 

Let us illustrate the concepts discussed in Sec. [TT] and 
in Sec. [Hi] with an established example, the Sherrington- 
Kirkpartick model.— Consider a system of N Ising spins 
<7j, i = 1,...,N, with phase space ilpf = {+1,-1}^, 
| Ojv | = 2 , subject to two-body interaction terms with cou- 
pling constants { Jy}i j=i n, J a — 0, Vi: 

E (W}) = -J2 J v ai<T j' (B1) 

where each Jy (= Jji) is independently Gaussian distributed 
with zero mean and standard deviation ~ 1 /N . Let us also as- 
sume that the dynamic processes are limited to single spin flip 
events, governed by Glauber dynamics. That is, a transition 
between two configurations C,C € Ojy differing by a single 
spin flip occurs with probability: 

e -0(E e ,-Ec)/2 

Pc ^ C ' = 2cosh[!3(Ec>-Ec)/2Y (B2) 
The SMF quantum Hamiltonian is then given by 4 

H = V 1 

SMF ^ i) 2co8b.\/3{Ec,-Ec)/2] 

X { e -KE e ,-Ec)/2 | C j ( C | _ | C ) ( C /| j j (B3) 

in the {\C}} basis, where (C,C) stands for pairs of configura- 
tions connected by a single spin flip, \C) = erf \C), 3i, 
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Using Eq. ( IB lb . and introducing the local fields hi 
^2 j Jij&ji we can write 

E e > - E e = 2 ^ ■l l ,n,n , = IhiGi 

3 

cosh[l3(E c , - Ec)/2] = cosh[/3/ii], 
and Eq. ( IB3b becomes 



Up to now, the variables cr, and hi are c-numbers. Introduc- 
ing the corresponding operators of and hi = J2j JijVj> one 
can simplify the notation and do away with the summation 
over C by observing that 



E 



i 



2 coshf/3/i. 



-Phial 



2 cosh[/3/i 4 



E 



i 



2 cosh[/3/i l 



|C) (C|af 



2 cosh[/3/i, 



(B4) 



-, (B5) 



where we used the fact that [Aj, of ] = (recall that J a = 0). 
We finally arrive at the expression, 



ffsMF = 



E 1 . {■ 



2 cosh[/3/i. 

Similarly, the GS wavefunction can be written as 4 

iv-o> = ^ f E e ^ c/2 i c ) 

10, 



<xf}. (B6) 



(B7) 



where Z = J2 C exp{^ £\ h i°i}- 

For convenience of notation, we will drop the : symbols 
from now onward, and all hi, of, and of will be understood 
as quantum mechanical operators. 

Given that the two operators e _/3/ti<T i and of in Eq. ( |B6| > 
do not commute, obtaining the spectrum of the system for all 
values of the parameter j3 is a tall order, in spite of knowing the 
exact GS wavefunction. In the following, we will investigate 
analytically the low and high temperature limits (f3 <C 1, and 
/3 » D- 



1. High-temperature behaviour 

The limit of /3 <C 1 can be studied by expanding the Hamil- 
tonian to first order in (3: 

h smf ~ ^ - 1 Yl w - \ E °f + (/3 2 )- ( B8 > 



Up to corrections of order /?, the GS is given by the (unique) 
eigenvector of of with eigenvalues +1 for all i, and the gap 
above it is equal to 1. According to Eq. (13. 2t and Eq. ( 13.3l l. 



this scenario corresponds to exponentially fast decaying auto- 
correlation functions, and vanishing Edwards-Anderson order 
parameter, as expected in the high-temperature phase of the 
dynamical classical system. 

[Note that the discussion so far applies in general to any 
classical Ising model with two-body interactions and Glauber 
dynamics.] 



2. Low-temperature behaviour 

The limit of (3 ^> 1 is more subtle, and it is convenient 
to approach it by looking at the GS wavefunction. We shall 
make use of working assumptions inspired by Parisi's work 
in Ref. [50f at low enough temperatures, the classical phase 
space effectively divides into "basins of influence" of the dif- 
ferent minima a = 1, . . . ,J\f of the free energy. That is, one 
can divide the partition function of the system Z = ^a, 
where all the relevant states for a given minimum a fall into 
the corresponding partial partition function Z a , and the ambi- 
guity in assigning all the states in between minima is immate- 
rial as their total weight in the partition function is negligible. 
For example, the reader could have in mind a 2D Ising fer- 
romagnet, where below T c one can safely divide the partition 
function into a positive and a negative magnetisation contri- 
bution (see for instance Sec. 3.1 in Ref. \^). 

Let us define the probability of being in a given basin, P a = 
Z a /Z {J2 a Pa = 1)- The GS wavefunction of the quantum 
SMF Hamiltonian can then be written as 



a 

|<M«^]Te-^/ 2 |C). 

V ^a, n ^ 



(B9) 



Notice that the partial wavefunctions \4> a ) have the same 
amplitudes as the GS wavefunction, but they involve only the 
relevant states in basin a. The assumption that the partition 
function can be divided into basins implies that local changes 
to a configuration C belonging to one basin will either produce 
a configuration in the same basin, or a configuration C whose 
weight is negligible, cxp[—/3Ec/2]/Z a <C e~ N . Since the 
off-diagonal term in the SMF Hamiltonian (with single spin- 
flip dynamics) is local, one can verify that this assumption re- 
sults in \<f) a ) being an approximate eigenvector of i?sMF with 
vanishing eigenvalue. 

The fact that no configuration state \C) appears in more 
than one Z a guarantees that the \(f> a ) are orthogonal to each 
other. They form therefore a basis of the low-lying manifold 
of eigenstates of the quantum system. 

We can then compute the Edwards-Anderson order param- 
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eter within this assumption [see Eq. ( I3.3b l : 

Af-1 2 

q E A(0)= |^0|O|^) 

n=l 

= E |(Vo|0hU| - \Wo\o\M 

n=0 

= ^o|o(e \4>n){4>n\\6\<k)- 

\n=0 j 

= Wo\ov&Wo)- {Mo Wo 



(V>o|0|Vo> 



(BIO) 



where the operator V is nothing but a projector onto the man- 
ifold of the Af low-lying states. We can therefore write it as 
V = Yl a =i \4>a) (0a |- Substituting into the previous equation, 
we arrive at the result 



A" 



9ea(C) « (Vol 6J2 10a) (0a| O |V>t>) " (Vol 6 |V>o) 



.V 



YPa\{<Pa\0\K 



N 



^p a (0 a |o|0 a ; 



where we used the fact that the operators O are necessarily 
diagonal in the configuration basis {\C}} (since they derive 
from classical observables O), and therefore (<p a \ O \ <f>b) = 0, 
Va ^ b. 

Notice that the expectation value (4> a \ O \<fi a ) corre- 
sponds to the ensemble average of the classical observ- 
able O in the ath basin of the free energy, (O) = 
J2 Cea O c e X p[-PE c }/Z a . Therefore, 



Af 



a=l 



A' 



(BID 



If we consider for example O = of, the second term vanishes 
in the absence of an explicit Z2 symmetry breaking term. By 
taking the average over all sites, 



M N , 2 ,2 A/" 

<?ea(0) ~ E Pq E AT = E Pfl ^ea, (B12) 



a— 1 i— 1 



we recover indeed the Edwards-Anderson order parameter in 
Parisi's formulation^! averaged over all basins. 



3. Griffiths singularities 

As noted in Sec. [TTTJ a finite value of <7ea implies that 
the local static susceptibility diverges in the low-temperature 
phase. However, this is a sufficient but not necessary condi- 
tion. For example, Griffiths singularities - characteristic of 
quantum disordered systems - can appear above the transi- 
tion, and cause the static susceptibility to diverge while ^ea 
remains zero^ 

Although investigating these issues is beyond the scope of 
the present paper, it is interesting to take a brief look at the 



specific example of the SK model, in the limit of j3 <C 1, 
Eq. (HI): 

iw ~ -~ E J *i°?°* - \ E < B13 > 



It is precisely the Hamiltonian of a quantum Ising spin glass 
in a transverse field, with mean field random interactions Jij. 
If it is the case that the small-/? approximation captures the 
physics of the full Hamiltonian i/sMF in the high-temperature 
phase j3 < j3 c , one would not expect any Griffiths singularities 
appearing in the disordered phase of this model. Indeed, in 
Eq. dB13t the local static susceptibility diverges only at the 
transition^ 

Although they may not play a role in the SK model, Grif- 
fiths singularities are likely to appear in other SMF Hamil- 
tonians associated with disordered dynamical systems (e.g., 
systems where Jij is short ranged)^ In such cases, it will 
be interesting to understand whether the Griffiths singulari- 
ties actually correspond to observable features in the original 
classical system, possibly relating to the presence of separate 
characteristic temperatures, as is the case for the dynamical, 
the thermodynamic one-step replica-symmetry-breaking and 
the full replica-symmetry-breaking transitions in the p > 3 
Ising spin glass model. 



Appendix C: Variational approach to the slowest relaxing modes 

Here we show how one can construct variationally the col- 
lapsing states (slowest relaxing modes) in the associated quan- 
tum system. We shall consider for simplicity the case of a 
classical nearest-neighbour Ising model, in which case the 
ground state in the quantum system is only two fold degener- 
ate below T c (see also Refs.Qj [l9|and references therein). The 
results we obtain are however more general, and in Sec. IVIBI 
we show how they can be used to find a subextensive set of 
collapsing state in an ideal glass system. 

The SMF Hamiltonian for the nearest-neighbour Ising 
model is a special example of the class of Hamiltonians dis- 
cussed in App.lBl namely where hi = J^2(ij) a j> J being the 
coupling constant in the classical system. Eq. ( IB6I ) reduces 
then to 



H 



SMF 



E 



2cosh[/3J£ fo)( 7j] 



(CI) 



and the relative GS wavefunction is 



Wo) 



\C) , (C2) 



where Z = ^ c exp{/3J^^ (J j (J i}- [Remark on the nota- 
tion: the sums over (ij) are over all nearest neighbors pairs of 
sites i, j, unless there is an outside sum over the site i, namely 
Si Y^(ij)> i R which case Y^Uj) stands for the sum over all js 
that are neighbors to the i site in the external sum J^i-l 
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Recall that, by construction, the GS [ipo) is annihilated by 
each term in curly brackets in Eq. dCll ). The existence of other 
state(s) asymptotically degenerate with \ipo) (in the ther- 
modynamic limit) means that 

(V'olV'i) = (C3a) 
ffsMF - A! (C3b) 
lim Ai = 0, (C3c) 

jV->oo 

where TV is the size of the system. 

While finding an exact expression for \tpi) is rather diffi- 
cult, one can attempt a variational approach. First of all, let us 
write as 

'«=?^^i>4 e>(c4) 

for a generic set of coefficients 0(C). The notation (. . .) th 
stands for a thermal average in the original classical system, 
namely 

<l°| 2 )th = |E |0(C)| 2 cxp{/3J^a^}. (C5) 

c (ij) I 



Notice that this thermal average is to be interpreted here as 
a mathematical summation over all spin configurations C, ir- 
respective, for instance, of the loss of ergodicity that occurs 
under spontaneous symmetry breaking (i.e., (af ) th = both 
above and below T c ). 

One can then show that (V'olV'i) is proportional to the ther- 
mal average (0) th . So long as we find a set of coefficients 
0(C) where (C) th — > for N — > oo, we satisfy the first con- 
dition in Eq. ( IC3b in the thermodynamic limit. That is, the 
condition of orthogonality in ( IC3b requires finding an observ- 
able O whose average over all configurations vanishes in the 
thermodynamic limit. (Notice that if (0) th = because of a 
symmetry, as is the case in the Ising model, then (il,<o\ipi) = 
even for finite TV.) 

Given (C) th = 0, one can immediately find an upper bound 
for the energy gap above the GS of the system, 

Ai<(V»i|5smf|^i>. (C6) 

Clearly the process can be iterated. Given a second ob- 
servable O' that satisfies (O')th = and (£>£>') th = 0, the 
state | -02) constructed from O' in the same way that was 
constructed from O is orthogonal to both \ipo) and \ipi), and 
another eigenstate must exist which differs in energy from the 
GSbyA 2 < (V> 2 |#smf|V2> 

Let us look in more detail at the structure of these upper 
bounds. Using the Hamiltonian in Eq. ( ICU . we obtain 



(0i | H SMF \^)= 1 £ O* (C) exp J J2 *K i 0(C) exp i ^ £ v'tf 
{M hhZ c,c> { 2 («> J [ 1 

= m^-z ? ? K )|a - °* (C) °N 2cLWs (M)0f] f (C8) 

1 i [|O(C)| 2 -O*(C)O(C fc )]exp{-/3JE (fc0 ^} f _ 

= (S^? ^hWE^f] cxp X J ^ at 

[\0(C)\ 2 -0*(C)0(C k )} oxp{-/3JE (fc ,)^^} ' 
2cosh[^JE<fci> °f] , 

E- « -> < C9 > 



A' 



/th 



where denotes the spin configuration that obtains from C 
upon changing the sign of af, . [Recall that, in presence of an 
overall summation over k, the sum Y^,lkl) a b° ve stands for a 
summation over all sites / that are nearest-neighbours to k.] 



Since the ratio between the exponential and the hyperbolic 
cosine inEq. (|C9t is of the formx(a;+l/a;) _1 with x positive, 



irrespective of the sign of ai, we can simplify the upper bound 

(\0(C)f-0*(C)0(C k )) th 



^E 

k 



(CIO) 



th 



Therefore, if we are interested in the lowest ly- 
ing eigenstates, we need to choose O such that 

E fc (|0(CI 2 -0*(C)0(C)) th «(|0| 2 ) th . 
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Specifically for systems that undergo a symmerty break- 
ing phase transition, such as the classical <i-dimensional Ising 
model with nearest-neighbour interactions (d > 1), a conve- 
nient choice for 0(C) that satisfies the condition (0) th = is 
an odd function of the order parameter, such as the magnetisa- 
tion of the system, M(C) = af (recall that (. . .) th stands 
for the mathematical ensemble average summed over all spin 
configurations C). 

In the case of the Ising model, a convenient choice to ob- 
tain a low energy bound is 0(C) = tanh[A/(C)]. Note 
that M(C) is the extensive magnetisation of the sample, and 
tanh 2 [M(C)] > ta nh 2 (l) so long as M (C) ^ 0. The denom- 
inator in Eq. dClOb is bounded from below by 



<|0| 2 ) th = (tanh 2 [M(C)]) 



th 



> tanh 2 (l) 

> tanh 2 (l), 



E 



{C:M(C)#0} 



-PE C 



-I3E C 



(Cll) 



where the ratio between the partition functions of the Ising 
model with and without the constraint M(C) ^ tends to 1 
in systems with an even number of sites, and it is identically 1 
in systems with an odd number of sites. 

We can thus focus only on the numerator in Eq. (IC 10b . 



|0(C)| 2 -0*(C)0(C fc ) 
= tanh 2 [M(C)] - 

= tanh 2 [M(C)) <! 



tanh[M(C)] tanh[M(C fe )] 
tanh[M(C fc )] 
tanh[M(C)] 



(C12) 



While in general this quantity can be negative, we know by 
construction that the GS energy of the SMF Hamiltonian is 
exactly zero, and therefore any upper bound to Ai must be 
non-negative. Under single spin flip dynamics, M(Cfc) = 
M(C) ± 2. Assuming M(C) > 2 and M(C k ) < M(C) for 
convenience, 

\0(C)\ 2 -0*(C)0(C k )< 
tanh[M(C fc )] 



~ 1 ~ tanh[M(C)] 

< 1 - tanh[M(C fc )] 

< 2e-' 2M(Ck \ 



(CI 3) 



Similarly for M(C) > 2 and M(C k ) > M(C), one obtains 

|0(C)| 2 - 0*{C)0(C k ) < 2e- 2M (°\ (C14) 

and analogous results hold for the two corresponding cases 
with M(C) < —2, with opposite sign in the exponent (i.e., 
e 2M(C) and e 2M(C fe ) ) respectively). 

Therefore, the thermal average appearing in the numerator 
of Eq. ( IClOb can be interpreted as the ratio between the par- 
tition functions of the classical Ising model with and without 
an additional energy term that penalises magnetised states, 



(\0(C)\ 2 -0*(C)0(C k )) 



th 



< 



E 



,-\M{C)\ p -/3E c 



J2 e-P E c 



While this ratio is finite in the paramagnetic phase, for 
T < T c it becomes exponentially small in the system size 



~ exp(— r yL d ~ 1 ), with 7 a non-universal constant. This ex- 
ponential penalty is the least one can incur [as opposed to 
a cost depending on the extensive magnetization and thus 
~ exp(— mL d )], and it is the cost of a domain wall that keeps 
M(C) ~ in the e _|M(c) term. So the bound on the gap Ai 
from this variational state is: 



Ai < 



0(L d ) T>T C 
L d e -iL d - 1 t<T c 



(C15) 



Notice that the bound for T > T c is rather loose, and we could 
have found a much better bound of 0(1) by other means. 
However, it is the bound for the ordered phase T < T c that 
is of interest here. The spectrum collapses exponentially in 
the system size for d > 1. 

The corresponding eigenvector for the Ising model is 



\i>i 



tanh[£\ of] ex P { IT E fo > °]at } 



c A /(tanh 2 E i( rf]) 



\C) 



Qualitatively, one can understand this result by noting that 
tanh[£\ cr,f]/y / (tanh 2 E i of ]) th - ±1, and the excited 
state I '0i ) is essentially the antisymmetric superposition of the 
positive and negative magnetisation valleys in the free energy. 
As expected, below T c the two valleys become ergodically 
disconnected, and it takes an exponential time in the size for 
the system to migrate from one valley to the other, leading to 
two distinct slow-relaxing modes, and therefore two lowest- 
lying states in the associated SMF Hamiltonian. For any finite 
size system, the actual GS is given by the symmetric (node- 
less) superposition of the two states, whereas the antisymmet- 
ric superposition lies at a slightly higher energy Ai. However, 
we expect (and indeed we just showed) that Ai — > in the 
thermodynamic limit N — > 00, for T < T c . 

We illustrated this constructive approach to a variational 
low-lying excited state in the simple case of a nearest neigh- 
bour Ising model. Given some a priori knowledge on the or- 
der parameter that distinguishes the free energy minima at low 
temperatures, it can be straightforwardly applied to find upper 
bounds to the lowest energy levels, say, in quantum Hamilto- 
nians derived from classical systems that undergo a dynamical 
transition. In the case of the Sherrington-Kirkpatrick model, 
briefly discussed in App. [B] this can be done only at a for- 
mal level, since we do not know explicitly the actual form 
of the non-local operators that select one valley over another. 
On the other hand, in Sec. [VI] we illustrate how the lowest 
energy states can be explicitly constructed, say, in kinetically 
constrained models without disorder. 

[The reader might be interested in comparing our deriva- 
tion of the excited states and of the upper bound to the spec- 
tral gap using a variational principle in quantum mechanics 
with the more elaborate but exact calculation of the relaxation 
spectrum using transition currents in Sec. 3. 2.1 in Ref.lH] 
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Appendix D: Quantum views on the gonihedric model 

In Sec.[vT]we introduced the generalised gonihedric model, 
and we discussed a few specific cases. Here we present two 
distinct dual descriptions of the isotropic limit, J, 



J, 



J zx = J, and we show how the relative SMF Hamiltonians 
relate to toric-code-like models in three dimensions (i.e., in- 
finitely massive Z2 quantum gauge theories). 



1. Gonihedric model 

Let us first recall the structure of the generalised gonihedric 
model. The Ising degrees of freedom <Xj = ±1 live on the sites 
of a cubic lattice, where we shall label the square plaquettes 
by p. The energy of the system can then be written as 



(Dl) 



p iGp 



where i 6 p label the four sites at the corners of p. Notice 
that flipping a plane of a spins does not change the energy 
in Eq. dDU . This symmetry of the model results in a mini- 
mal degeneracy of isoenergetic configurations that scales as 

23L-2 

Assuming Glauber single spin flip dynamics, a few alge- 
braic steps 4 - as in the example discussed in Sec.|Vl]lead to the 
associated SMF Hamiltonian 



H, 



(o) 

SMF 



E 



1 



2 cosh 

AE i= 2j ]r 









| cxp 


2 





(D2) 



{p:i£p} j&P 



2. Gonihedric model in the bond-dual spin language 

An alternative description of the same classical system can 
be formulated in terms of Ising spins living on the bonds 6 of 
the lattice, 9^ = 0j_(&)<Xj + n,), where i±(b) are the two sites 
adjacent to bond b, as illustrated in Fig. [7] In Sec. I VI A~T1 dis- 




FIG. 7. Illustration of the mapping from a to 9 spins. 

cussing the anisotropic gonihedric model, we used only the 
9 spins on the vertical bonds of the cubic lattice. Indeed, 



with the addition of a horizontal plane of a spins, they are 
sufficient to describe the entire system. Here we define in- 
stead a 9 spin on each bond of the lattice. Whilst the nota- 
tion becomes highly redundant, it demonstrates an interesting 
analogy between the SMF Hamiltonian for the isotropic go- 
nihedric model and a class of topologically ordered systems 
called toric code models.— 

In this language, given a plaquette with corners i\, . . . , 14, 
and bonds 61 , . . . , 64, the corresponding 4-body energy contri- 
bution cr, 1 (jj 2 (Tj 3 o'j 4 becomes a 2-body term that can be writ- 
ten as 9^9^ or 9^9^. However, not all {9^} configurations 
are allowed by the mapping. Indeed, for every plaquette we 
have the gauge constraint 9i >1 9j, 2 9f, 3 0j, 4 = 1. In addition, the 
product of all the 9 spins along a winding line on the lattice 
(periodic boundary conditions are assumed) must be positive. 
Even if there are 3N 9 spins for a system of TV a spins (i.e., N 
sites of the cubic lattice), the constraints reduce the number of 
independent 6s, and ultimately one can see that the mapping 
is l-to-2, since one can obtain the value of any a spin given 
all the 9s plus one reference a. The energy of the system can 
be written as 



2 ^ 

[[66']] 



'bVV 



(D3) 



where the notation [[66']] stands for two parallel bonds 6 and 
6' belonging to the same plaquette. 

Note that Glauber single spin flip dynamics in the a spins 
translates into a 'star' flipping operation, that changes the sign 
of all 6 9t spins adjacent to a common site i (6 £ i). The 
associated SMF Hamiltonian assumes the form 



H, 



(i) 

SMF 



E 



1 



2 cosh 



exp 



PAEi 



n< 

6Gi 



AE > = J E E «■ 

6Gi [[66']] 



(D4) 



Notice that all terms in the Hamiltonian commute with the 
constraints, namely the product of all 9 Z components of the 
spins around a plaquette p is the identity, Ylbep ®l = 1> an< ^ 
that the product of 9 Z along any winding loop 7 also gives 
the identity, Yibey ®b = We can therefore extend the 
zero-temperature SMF quantum system to the unconstrained 
Hilbert space generated by the spin- 1/2 9 degrees of freedom, 
provided that we add an appropriately large energy cost to 
states that violate the constraints: 



H, 



(i) 

SMF 



E 



1 



2 cosh 



En< 

p 6Sp 



n 

6e7i 



exp 



n< 

b£i 



n« 



n 

6G7* 



6 • 



(D5) 



where j x , j y and j z are three generic winding loops along the 
x, y and z direction, respectively. 

Is is interesting to notice that the limit J = in Eq. dP5t is 
nothing but Kitaev's toric code Hamiltonian in 3D ) 16 - 54 ' 55 with 
the addition of explicit non-local operators (the products along 
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the 7 Q winding loops) that select one out of the 8 topological 
sectors. The (gapped) topological phase corresponds to the 
high-temperature limit of the classical system, and it survives 
for small values of j3 up to a phase transition (the transition in 
the classical gonihedric model 38 ) where topological order is 
lost, and yet neither of the two neighbouring phases appears 
to be captured by a local order parameter. 

3. Gonihedric model in the face-dual spin language 

Finally, there is another useful description of the isotropic 
gonihedric model, formulated in terms of Ising spins living on 
the plaquettes p of the lattice, S p = Yiiep as illustrated in 
Fig. [8] Note that the S spins live on the bonds of the cubic 




FIG. 8. Illustration of the mapping from a to S spins. The S spins 
living at the plaquette centres in the original lattice can be also seen 
as bond degrees of freedom in the dual cubic lattice (thick light lines 
in the figure) formed by the centres s of the unit cubic cells in the 
original lattice. The four S spins belonging to a planar star on the 
dual lattice, centred at s and perpendicular to the x direction, are 
highlighted by the dashed blue lines. 

lattice formed by the centres of the unit cubic cells in the orig- 
inal lattice. In this language, the energy becomes a trivial sum 
of 1-body terms 

E = -JJ2S P , (D6) 

p 

As above, not all {S p } configurations are allowed by the 
mapping. For every unit cubic cell in the lattice there are three 
constraints that must be satisfied. Namely, the product of the 
four plaquette spins on the faces of the cell parallel to a lat- 
tice direction must be equal to the identity, for each direction 
x, y, and z, as illustrated in Fig. [9] (see Ref. |56| for a re- 
lated duality transformation). In the dual cubic lattice, these 
identities correspond to the condition Il{pes P j_ a } S p = +1, 
where s denotes a generic dual lattice site, p E s labels the 6 
plaquette sites on the adjacent bonds, and p _L a means that 
the dual bond on which p lives is perpendicular to the lattice 
direction a = x,y, z. In addition, the product of all plaquette 
spins S p along any strip (equivalent to a spin ladder, in the 
original a spins) winding around the system is also equal to 
1. In the end, one can see that the mapping is l-to-(3i — 2), 




FIG. 9. Illustration of the local constraints in the dual S p spin lan- 
guange in the y direction: the product of the S p spins on the shaded 
plaquettes must always be equal to 1. Similarly for the x and z di- 
rections. 

since we can obtain the value of any a spin given all the S 
spins plus all the a spins along three reference straight lines 
winding around the system, and crossing all at one point in 
the lattice. Note that 3L — 2 is precisely the degeneracy of the 
lowest energy states in the gonihedric model. Indeed, chang- 
ing the sign of all the spins belonging to the same lattice plane 
is an exact symmetry of the classical energy, and it allows 
one to construct 3L — 2 distinct lowest energy configurations 
starting from, say, the fully magnetised one {er.; = +1}. All 
such configurations correspond to the same fully magnetised 
{S p = +1} plaquette spin configuration, in agreement with 
the l-to-(3L — 2) nature of the mapping. 

The Glauber single spin flip dynamics in the a spins trans- 
lates into a flipping operation that changes the sign of all the 
Sp spins on plaquettes that have the flipped spin er.; at one of 
their corners ({p : i £ p}). Alternatively, this can be seen as 
the operation that flips all S p living on the 12 edges of a cubic 
unit cell in the dual lattice (thick light green lines in Fig. [8]), 
centred at i. The associated SMF Hamiltonian assumes thus 
the form 

*-?5^H-^];n?} 

AE i = 2J J2 S p- ( D7 > 

{P ■■ «Gp} 

All terms in the Hamiltonian commute with the constraints, 
namely the products of the S z components of the spins around 
planar stars and straight winding lines 7 (perpendicular to the 
centres of plaquettes in the original lattice) being the iden- 
tity. We can therefore extend the zero-temperature SMF quan- 
tum system to the unconstrained Hilbert space generated by 
the spin- 1/2 S degrees of freedom, provided that we add an 
appropriately large energy cost to states that violate the con- 
straints: 




E E II ^-EIK*> < D8 > 

a=x,y,z s {p£s:p±a} 7 P£7 
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where s labels the sites in the dual cubic lattice. 

In the limit of J = 0, Eq. dD8l ) reduces to a Z 2 lattice gauge 
theory different from the one obtained in the previous section. 
This gauge theory has a peculiar topological degeneracy that 
scales with the size of the system to the power 2/3 (namely, 
the sectors are identified by the eigenvalues of the ~ 3L 2 7 
winding loop operators). As in the case of the 6 spin mapping, 
Eq. ( ID8b contains non-local terms that select a unique topo- 
logical sector. The (gapped) topological phase corresponds to 
the high-temperature limit of the classical system, and it sur- 
vives for small values of j3. 



Now, the /?->oo limit can be recognized as the trivial fully 
magnetised state with S p = +1 for all p, which is unique 
and non-degenerate (in the S p variables of this face-dual de- 
scription). We can distinguish between the two phases using 
the topological entropy defined in Refs. [261271 Here we find 
Sto P o(/3 = 0) ^ and Sto P o(P -> 00) = 0. We can there- 
fore prove the existence of a phase transition separating these 
two phases at some critical temperature j3 c , in support of the 
current evidence based on numerical simulations and cluster 
mean field arguments^ (See also Ref. H3 for an alternative 
proof of the transition.) 
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